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We use the recently extended superposition T-matrix method to calculate scattering and absorption properties of 
micrometer-sized water droplets contaminated by black carbon. Our numerically exact results reveal that, depend- 
ing on the mode of soot-water mixing, the soot specific absorption can vary by a factor exceeding 6.5. The specific 
absorption is maximized when the soot material is quasi-uniformly distributed throughout the droplet interior in the 
form of numerous small monomers. The range of mixing scenarios captured by our computations implies a wide 
range of remote sensing and radiation budget implications of the presence of black carbon in liquid-water 
clouds. We show that the popular Maxwell-Garnett effective-medium approximation can be used to calculate 
the optical cross sections, single-scattering albedo, and asymmetry parameter for the quasi-uniform mixing scenario, 
but is likely to fail in application to other mixing scenarios and in computations of the elements of the scattering 
matrix. © 2014 Optical Society of America 
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It is widely recognized that contamination by black 
carbon may cause substantial changes in the optical 
and radiative properties of liquid-water clouds, with 
strong implications for Earth’s climate system (see, 
e.g., Refs. [1-3] and references therein). However, accu- 
rate quantification of these changes had been a very 
challenging problem since direct computer solutions of 
the Maxwell equations for realistic models of microm- 
eter-sized droplets with various types of soot impurities 
were impracticable [4,5]. Instead, the computations of 
electromagnetic scattering by such particles had been 
based on the simplified model of a droplet with only 
one inclusion [1] or on approximate approaches, such 
as various effective-medium theories (EMTs) [2,3,6] 
and the Monte Carlo ray-tracing technique [7,8]. This un- 
fortunate situation has changed owing to the recent ex- 
tension of the superposition T-matrix method (STMM) 
[9,10] to arbitrarily clustered and nested spherical do- 
mains [11]. Hence, the principal objective of this Letter 
is to perform the first quantitative comparison of scatter- 
ing and absorption properties of different types of drop- 
let-soot mixtures. We will also analyze to what extent 
these properties can be reproduced by a popular EMT, 
viz., the Maxwell-Garnett approximation (MGA) [6]. 

To cover a wide range of potential mixing scenarios, 
we consider the following 11 models (sketched in Fig. E) 
while keeping the wavelength fixed at A = 0.55 pm. 

Model A. A fluffy soot aggregate composed of N = 75 
identical spherical monomers is internally mixed with 
(i.e., is completely imbedded in [12]) a spherical water 
droplet such that the distance between the droplet center 
and the geometrical center of the fractal is equal to d. The 
droplet radius, R, is such that the droplet size parameter 
X = 2nR/X is equal to 40, while the droplet refractive 
index is 1.33. The morphology of the soot aggregate is 
parameterized by the standard statistical scaling law 
N = ko(R n /r) D f [13], where r is the monomer radius, 


k 0 = 1.19 is the fractal prefactor, D f = 1.82 is the fractal 
dimension, and R g , called the radius of gyration, is a mea- 
sure of the overall aggregate radius. The black carbon 
refractive index is fixed at 1.95 + 70.79, while the mono- 
mer radius is fixed at the value corresponding to the 
monomer size parameters = 2 nr /X = 1 (note that these 
numbers do not necessarily represent the most typical 
experimental or modeling values; cf. [14-16]). All scatter- 
ing and absorption properties are averaged over the 
uniform orientation distribution of the resulting inhomo- 
geneous particle. 

Model B. As in Model A, but the fluffy soot aggregate 
collapses into a compact cluster with D f = 3. 
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Fig. 1. Models of soot-contaminated water droplets. 
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Model C. As in Model B, but the compact soot cluster is 
replaced by a large homogeneous spherical soot inclu- 
sion having the same volume as the soot cluster. 

Model D. As in Model A, but the spherical water droplet 
and the fluffy soot aggregate are externally mixed [12], 
i.e., are widely separated and act as “independent scat- 
tered”. The radius of the droplet is adjusted to keep 
the combined volume of the water and black carbon 
fixed, the corresponding droplet size parameter being 
39.984. All scattering and absorption properties of the 
soot cluster are averaged over the uniform orientation 
distribution. 

Model E. As in Model B, but the spherical water droplet 
and the compact soot aggregate are externally mixed and 
the radius of the droplet is adjusted. All scattering 
and absorption properties of the soot cluster are aver- 
aged over the uniform orientation distribution. 

Model F. As in Model C, but the spherical water droplet 
and the homogeneous spherical soot particle are exter- 
nally mixed and the radius of the droplet is adjusted. 

Model G. As in Model D, but the fluffy soot aggregate 
and the droplet are semi-extemally mixed [12], i.e., the 
aggregate is attached externally to the surface of the 
droplet. 

Model H. As in Model E, but the compact soot aggre- 
gate and the droplet are semi-extemally mixed. 

Model I. As in Model F, but the homogeneous spherical 
soot particle and the droplet are semi-extemally mixed. 

Model J. As in Model A, but the 75 soot monomers 
are quasi-randomly distributed throughout the droplet 
interior rather than form a cluster. 

Model K. A homogeneous droplet has a size parameter 
of 40, and its refractive index is computed according to 
the MGA for the same volume of soot material as in 
Models A-J. 

Figure 2 depicts the orientation-averaged absorption 
cross section, C abs , for Models A-C as a function of 
the size-parameter distance, 2nd IX. It is obvious that, 
for the same total amount of black carbon, the cumula- 
tive absorption (i) strongly depends on the specific mor- 
phology of the soot inclusion and (ii) is maximal when 
the soot inclusion is located sufficiently far from the 
droplet boundary. In geometrical optics terms, distancing 
the inclusion from the boundary and increasing its geo- 
metric cross section serves to maximize the probability 
of the inclusion to “intercept” refracted and internally 
reflected “rays”. This apparently explains both traits 
visible in Fig. 2, at least qualitatively. 

The first trait is further illustrated by the third column 
of Table 1, which shows the absorption cross section 
for Models A-C averaged over all relevant distances, 
d. This column also demonstrates that absorption be- 
comes minimal when the entire soot volume is concen- 
trated in one homogeneous spherical particle that is 
mixed with the droplet externally or semi-extemally. On 
the other hand, breaking a soot fractal into individual 
monomers and distributing them quasi-uniformly 
throughout the droplet yields the greatest C abs value. It 
is quite remarkable that keeping the amount of soot fixed 
but switching from Model F to Model J serves to increase 
the specific soot absorption by a factor exceeding 6.5. 

Table 1 also lists the respective values of the extinction 
cross section, C ext , scattering cross section, C sca , 



Fig. 2. Absorption cross section versus 2nd/ X for Models A-C. 

asymmetry parameter, g , and single-scattering albedo, 
tu. It is seen that these integral characteristics are rela- 
tively less affected by the soot morphology and type of 
mixing. However, the single-scattering albedo difference 
between Models F (or I) and J is quite significant since it 
can result in a substantial change of the reflectivity of an 


Table 1. Integral Radiative Characteristics of Water 
Droplets Mixed with Black Carbon 


Model 

Cext (l-™ 2 ) 

Cabs (W) 

Csca Om 2 ) 

g 

TU 

A 

77.9846 

2.4649 

75.5197 

0.8395 

0.9684 

B 

77.3582 

1.3594 

75.9988 

0.8398 

0.9824 

C 

77.0620 

0.9443 

76.1177 

0.8396 

0.9877 

D 

79.7094 

1.5891 

78.1203 

0.8362 

0.9801 

E 

78.4154 

0.8850 

77.5304 

0.8376 

0.9887 

F 

77.8879 

0.5753 

77.3126 

0.8370 

0.9926 

G 

78.6325 

1.5835 

77.0489 

0.8373 

0.9799 

H 

77.4591 

0.8757 

76.5834 

0.8387 

0.9887 

I 

77.1285 

0.5763 

76.5522 

0.8380 

0.9925 

J 

77.0129 

3.7469 

73.2660 

0.8412 

0.9513 

K 

76.6184 

3.9245 

72.6939 

0.8510 

0.9488 
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optically thick liquid-water cloud. Indeed, the cloud al- 
bedo is largely defined by the so-called single-scattering 
co-albedo, 1 - m, rather than by m. Table XIX on page 340 
of [17] shows that switching from Model F to model J can 
be expected to reduce the spherical albedo of an opti- 
cally semi-infinite cloud by almost a factor of 2. Another 
important result is that the pairs of external and semi- 
external mixing scenarios (D,G), (E,H), and (F,I) yield 
nearly identical results. 

The last row of Table 1 shows the result of applying the 
MGA. It is quite obvious that the MGA values of the 
absorption cross section and single-scattering albedo 
reproduce reasonably well the STMM results for 
Model J but not those for all the other models. 

The angular distribution and the polarization state of 
the scattered light can be conveniently described in 
terms of the normalized Stokes scattering matrix, 


F(0) = 
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where 0 e [0°, 180°] is the angle between the incidence 
and scattering directions [18]. The (1,1) element of 
F(0) is the conventional phase function normalized 
according to 


d0 sin 0 ^(0) = 1. (2) 

Two unique properties of a homogeneous, perfectly 
spherical particle are the identities 

a 2 (0)/ai(0) = l and a 3 (0)/a 4 (0) = 1. (3) 

Any EMT always reproduces these identities by defini- 
tion, provided that the surface of the host is perfectly 
spherical. However, Fig. 3(a) reveals a fundamental vio- 
lation of the first identity by the numerically exact STMM 
results and thereby illustrates an essential limitedness of 
the concept of effective refractive index. The violation of 
the second identity (not shown) is equally significant. 
Interestingly, the angular locations of all minima in 
Fig. 3(a) appear to be the same irrespective of the mixing 
scenario. 

Figures 3(b) and 3(c) show STMM results for anl = 
20 droplet quasi-uniformly filled with N soot inclusions, 
each having a size parameter x = 1. Fig. 3(b) shows that, 
with increasing N, the deviation of the ratio a 2 (0) /<Xi(0) 
from unity first increases, then saturates, and finally 
starts to decrease. This behavior can be explained quali- 
tatively by competing effects of increasing “multiple scat- 
tering” and increasing absorption. The former serves to 
reduce the ratio whereas the latter serves 

to suppress the contribution of scattered light emerging 
from the droplet interior and thereby enhances the 
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orientations and, for Models A and C, over the applicable ranges of d. (b), (c) STMM and MGA results for a water droplet filled quasi- 
uniformly with N small black carbon inclusions, (d), (e) MGA errors as functions of the droplet size parameter or the number of 
black carbon inclusions. 
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relative contribution of the external reflection satisfying 
the first identity of Eq. (3). In general, Figs. 3(a) and 3(b) 
demonstrate the potential of using measurements of the 
ratio a 2 (0) /ai(&) as an optical characterization tool 
highly sensitive to nonspherical and/or inhomogeneous 
particle morphology (cf. [19,20]). 

The phase function results plotted in Fig. 3(c) for the 
same X = 20 droplet show that increasing N serves 
to make a^©) progressively featureless at side- and 
back-scattering angles and to suppress the backscatter- 
ing phase function values. The effect of applying the 
MGA is qualitatively similar, but results in quantitative 
errors exceeding a factor of 3 in the case of the droplet 
with 800 soot inclusions. 

Figures 3(d) and 3(e) further quantify MGA errors in 
integral radiometric characteristics computed for Model 
J assuming a size parameter of the soot inclusions fixed 
at x = 1. In Fig. 3(d) , the droplet size parameter, X , is 
increased from 10 to 46.42, while the number of inclu- 
sions is kept fixed at 100, thereby implying a decrease 
of the black carbon volume fraction from 10% to 0.1%. 
In Fig. 3(e) , the droplet size parameter X is fixed at 
20, while the number of inclusions, N, grows from 5 to 
800. The black carbon volume fraction thus increases 
from essentially zero to 10%. It is seen that, in all cases, 
the MGA errors do not exceed ~10%. 

In summary, our numerically exact STMM results 
reveal a very strong dependence of the soot specific 
absorption on the mode of soot-water mixing in 
liquid-water clouds contaminated by black carbon. The 
absorption is maximized when the soot material gets 
quasi-uniformly distributed throughout the droplet 
interior in the form of numerous small monomers (Model 
J). The absorption is minimal when black carbon exists in 
the form of larger homogeneous particles mixed with 
cloud droplets externally or semi-extemally (Models F 
and I). Similar conclusions were reached in a recent 
modeling study of mixtures of soot aggregates with small 
sulfate aerosols based on the discrete dipole approxima- 
tion [21]. The presence of soot has other noticeable man- 
ifestations, such as its effects on the single-scattering 
co-albedo and the elements of the scattering matrix. It is 
beyond the scope of this Letter to discuss the relative 
plausibility of the mixing scenarios represented by 
Models A-J. However, the morphological range captured 
by these scenarios implies a wide range of remote sens- 
ing and radiation budget implications of the presence of 
black carbon in liquid-water clouds. 

An important by-product of our study is the first 
quantitative analysis of the accuracy of the MGA as 
applied to micrometer-sized water droplets mixed with 
black carbon. It appears that the only worthwhile 
application of the MGA is the calculation of the optical 
cross sections, single-scattering albedo, and asymmetry 


parameter for the quasi-uniform mixing scenario des- 
cribed by Model J. 

This material is based upon work supported by the 
NASA Remote Sensing Theory Program managed by 
Lucia Tsaoussi and the NASA Radiation Sciences Pro- 
gram managed by Hal Maring. The majority of numerical 
results were obtained with the “Discover” supercom- 
puter at the NASA Center for Climate Simulation. 
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